Supersolidity in quantum films adsorbed on graphene and grapliite 
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Using quantum Monte Carlo we have studied the superfluid density of the first layer of *H[e and 
H2 adsorbed on graphene and graphite. Our main focus has been on the equilibrium ground state of 
the system, which corresponds to a registered y/3 x y/S phase. The perfect solid phase of H2 shows 
no superfluid signal whereas *He has a finite but small superfluid fraction (0.67%). The introduction 
of vacancies in the crystal makes the superfluidity increase, showing values as large as 14 % in '*He 
without destroying the spatial solid order. 
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Supersolid state of matter is a fascinating possibility 
that has long attracted interest from both theoretical and 
experimental viewpoints. ^'^ The simultaneous existence 
of spatial lattice order and off-diagonal long range or- 
der defining a supersolid is rather counterintuitive and 
only a theoretical entelechy up to recent times. The old 
theoretical ideas put forward by the pioneering works of 
Andreev and Lifshitz-^. and Chester— and Reatto^ have re- 
vived dramatically since the experimental findings of Kim 
and Chan^ in 2004 on the evidence of non-conventional 
moment of inertia in solid ^He below some characteristic 
temperature around 100 mK. Whether these and some 
other similar experiments carried out by other teams are 
an unambiguous proof of the existence of supersolidity or 
not is still a matter of debate. 

At present, much less is known on possible super- 
solid scenarios in two-dimensional (2D)1 or quasi-two- 
dimensional (Q2D) solid *He. Helium atoms, when ad- 
sorbed on graphite or graphene, arrange sequentially in 
stacking layers that can be considered as nearly two- 
dimensional systems. * Interestingly, recent experiments 
carried out by Nyeki et al.^ in the second layer of ^He ad- 
sorbed on graphite point to the existence of a (a/7 x \/7) 
commensurate solid phase that exhibits superfluid frac- 
tions as large as 20%. However, recent path integral 
Monte Carlo (PIMC) work^ does not seem to support 
the existence of this phase which was first predicted in 
simulations where the first layer of ^He atoms was con- 
sidered as inert Nevertheless, the possibility of having 
a supersolid in a 2D environment, supported by Nyeki et 
al.'s observations)^ opens new and exciting avenues for 
the analysis of this new state of matter. 

Recently, we have calculated the zero-temperature 
phase diagram of the first layer of ^He^ and H^^^ ad- 
sorbed on graphene and graphite using quantum Monte 
Carlo (QMC) methods. Our results predict that the 
equilibrium ground state of both systems is a 
commensurate phase, conclusion which, in the case of 
graphite, is in agreement with other low-temperature 
simulations and experimental dataJ^^— The aim of our 
previous studies was basically the determination of the 



energies of the different possible phases to draw the phase 
diagram and not the study of off-diagonal long range or- 
der and/or superfluidity. To this end, we used a non- 
symmetric wave function for describing the solid phases, 
an approach which obviously hinders any insight on the 
properties directly related to their Bose-Einstein statis- 
tics but which guarantees accurate evaluation of the en- 
ergies due to the low-rate interparticle exchange. In this 
work, we are mainly concerned with the possibility of 
supersolidity in the first layer of ''He and H2 adsorbed 
on graphene and graphite so that our methodology has 
been changed accordingly. Besides the characterization 
of the equilibrium ground-state phases, we have also an- 
alyzed the influence of vacancies in the superfluid re- 
sponse and energy of these films since point defects are 
indeed observed during quantum layer nucleation on car- 
bon surfaces. It is worth mentioning that previous low- 
temperature attempts using PIMC have not found any 
signal of superfluidity in these layers^°'^^ although possi- 
ble supersolidity induced by defects was not analyzed. 

Since we are interested in possible ground-state super- 
solid phases of the first ''He and H2 layers adsorbed on 
graphene and graphite, we use the diffusion Monte Carlo 
(DMC) method that, working at zero temperature, solves 
stochastically the Schrodinger equation of the A^-particle 
system in an essentially exact way (within some statis- 
tical uncertainties)!^ Zero-temperature approaches are 
specially adequate for these systems since the onset tem- 
peratures for supersolidity are, at least in bulk solids, 
very small (50-100 mK) thus making finite-temperature 
approaches like PIMC of difficult use. The Hamilto- 
nian of the system includes both interparticle {V(r)) and 
particle-substrate {U{r)) interactions. 
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with capital indexes running on carbon substrate atoms 
and normal indexes on He atoms or H2 molecules. The 
''He-''He and H2-H2 interactions are modeled by the ac- 
curate Aziz II^'^ and Silvera-Goldmani^i potentials, re- 
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FIG. 1. (Color online) (a) Energy per *He atom with differ- 
ent number of vacancies, plotted as a function of the particle 
density; (b) Same as in (a) but for H2. In both cases the 
energy per particle of the commensurate phase with no va- 
cancies {E/N)o is subtracted. Solid line and squares (dotted 
line and diamonds) stand for energies when the adsorbate is 
graphene (graphite). 



spectively. The adsorbate surface potential is obtained 
by summing all the pair atom(molecule)-adsorbate inter- 
actions modeled by Lennard- Jones potentials. By sum- 
ming up all the latter pair interactions, we introduce in 
the description of the system the necessary corrugation 
to observe that a commensurate solid phase is effectively 
preferred for being the ground stateJ^ii^ 

As it is usual in DMC, we introduce a trial wave 
function for importance sampling which improves the 
variance of the statistical estimation. Our variational 
model contains basic information: it is zero at shorter 
distances than the (hard-)core of the potential and be- 
comes constant at large distances (Jastrow wave func- 
tion); it is symmetric under the exchange of particles 
(Bose-Einstein statistics), and localizes particles in pre- 
ferred points (sites) for solid phases. Explicitly)^ 
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In Eq. (121), R = {ri, . . . , rftr}, /(r) is a two-body Jastrow 
correlation factor, g(rji) = exp[— a(ri — r/)^], and Ncr is 



the number of lattice sites of the selected crystal struc- 
ture. Model wave function ([2]) fulfills simultaneously spa- 
tial solid order and exchange-particle symmetry avoiding 
the numerically unworkable use of permanents on top of 
the non-symmetric Nosanow- Jastrow wave function. The 
value of the variational parameters in are the same as 
the ones reported in Refs. [H andlTsl. 

We have focused our attention in the ^/Sx \/3 commen- 
surate phase of ^He and H2 since this is the equilibrium 
ground state structure and best candidate for exhibiting 
supersolid behavior due to its quite low density. The sim- 
ulation cell is a rectangle of fixed dimensions 44.27 x 42.60 
where the number of sites of the registered phase is 
A'cr = 120 and the number of particles is 115 < N < 120 
,i.e., we consider up to a maximum of five vacancies. The 
dependence of the energy per particle with the number of 
vacancies, reported as a function of the particle density, is 
shown in Fig. 1 for ^He and H2 on top of graphene and 
graphite. For comparison purposes, those curves have 
been displaced downwards an energy shift {E/N)q which 
corresponds to the energy of the perfect structures re- 
ported in Refs. [l^ and[l3l For the present analysis it is 
enough to mention that the binding energy of H2 is ap- 
proximately 2.5 times larger than that of ^He and that 
in both systems the interactions with graphite are about 
10 % more attractive than with graphene. 

As shown in Fig 1, the binding energies per particle 
increase systematically with the number of vacancies so 
the equilibrium ground state always corresponds to the 
perfect commensurate lattice. The variation of the acti- 
vation vacancy energy is linear in both cases and slightly 
larger in graphite where the adsorption energies are in- 
variably ~ 10 % larger than in graphene. Comparing the 
results obtained for ''He and H2, one can see that the en- 
ergy difference E/N — {E/N)o is roughly 10 times larger 
for hydrogen however, as the absolute values for both 
solids are so different, it is better to make a comparison 
in relative terms: the energy cost of creating 5 vacancies 
in "^He is 0.15 % of the ground-state energy while in H2 
is 0.43 %. When the number of vacancies increases, and 
therefore the particle density decreases, the equation of 
state of the solid with vacancies approaches the equation 
of state of the metastable liquid phase. This analysis is 
specially interesting when graphene is the substrate since 
the difference in binding energy at the equilibrium point 
of the liquid and the one at the commensurate solid phase 
is very tiny, nearly 4 times smaller than in graphite. For 
both ''He and H2, all the energies shown in Fig. 1 are 
below the ones calculated for the corresponding liquid 
phases at the same densities. If the linear behavior ob- 
served in the figure is extrapolated to smaller densities, 
one can see that in ^He the intersection with the liquid 
equation of state is produced at density CTc = 0.058 A^^, 
which would correspond to 10 vacancies in the simulation 
box. In this case, CTc is larger than the equilibrium density 
of the liquid ctq = 0.044 A^^ thus the crossing point will 
appear at finite pressure. In contrast, the intersection 
with the liquid equation of state in hydrogen appears at 
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FIG. 2. (Color online) DMC estimation of the superfluid den- 
sity. Solid, dotted, and dashed lines stand for commensurate 
*He, incommensurate *He, and commensurate H2 phases on 
top of graphene, respectively. 



density Uc = 0.052 A~^, which is smaller than the equi- 
librium point (Jo = 0.059 A^^ and thus corresponds to 
the metastable negative pressure regime. 

A discussion on possible supersolid phases of helium 
and hydrogen on top of graphene or graphite requires the 
estimation of the superfluid density fraction of the differ- 
ent solid phases that appear in their respective phase dia- 
grams. In quantum Monte Carlo, the superfluid fraction 
can be computed by sampling the winding number— 
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with T the imaginary time and /3 = T^^. In the limit of 
zero temperature /3 — ^ cxi, the winding number ([3]) turns 
to the diffusion coefficient of the center of mass of the 
N particles (Rcm) in the simulation box with periodic 
boundary conditions;^ 
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where a — N/2dDQ, with d the number of dimensions 
(d = 2 in the present case), Dq — h^/2m, and Ds{t) — 
{{RcMir) - HcA/(0))2). The diffusion coefficient D,{t) 
can be calculated using the DMC method, in which the 
imaginary time evolves in a continuous way, and it can 
be proved that in the asymptotic regime this estimator 
is unbiased with respect to the particular choice of the 
trial wave function used for importance sampling (within 
a specified physical phase). According to the typical sta- 
tistical noise in DMC simulations, the resolution of this 
estimator is much lower (^ 1 x 10^^) than in PIMC, 
where superfiuid signals below ^ 1 x 10^^ arc extremely 
difficult to be measured. 

In Fig. 2, results for aDs{T) as a function of the imagi- 
nary time T are shown. As obvious from its definition (jj]). 
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FIG. 3. (Color online) (a) Superfluid fraction of the *He com- 
mensurate solid phase with different number of vacancies ; 
(b) Same as in (a) but for H2. Solid line and squares (dotted 
line and diamonds) stand for superfluid fractions when the 
adsorbate is graphene (graphite). 



a finite superfluid fraction appears as a finite slope in the 
long-time behavior of the diffusion coefficient Dg (r) , their 
particular values being not relevant and fairly dependent 
of the kind of system under study. Our results for the 
perfect (no vacancies) solid phases plotted in Fig. 2 show 
different behaviors depending on the system and solid 
phases considered. Diffusion coefficient results for long 
time T obtained for the x ^/S commensurate phase 
of ''He on top of graphene (and also in graphite) show a 
small but clear slope {ps/p = 0.0067(1)) that contrasts 
with the null r-variation observed in a incommensurate 
phase of density 0.0999 A~^. Interestingly, simulations 
performed in the commensurate phase of H2 indicate zero 
superfluid fraction thus areal density must be ruled out 
as the only cause behind supersolidity in quantum films. 
Moreover, as it was shown in the DMC calculation of 2D 
and quasi-2D "^He in Ref. [2^, the superfiuid fraction of 
a purely two-dimensional crystal is zero even at densi- 
ties below the a/3 x a/3 phase, and finite superfiuidity 
emerges only with the opening of a transverse direction 
that particles can explore. According to these previous 
results, the zero signal observed in H2 can be explained 
in terms of transverse motion frustration resulting from 
intense molecular binding to carbon surfaces. 



4 




4 8 12 16 

x[A] 




FIG. 4. (Color online) Top: Snapshot of the probability den- 
sity (crosses) of commensurate solid *He on top of graphene. 
Bottom: Same as in the top panel but with five less particles 
in the simulation box. In both figures, squares are the sites 
of the perfect crystal. 



One of the aims of this work has been the study of 
the influence of the number of vacancies in the superfluid 
fraction of quantum sohd layers. To this end, we have 
calculated Ps/ p for the commensurate phase of ^He and 
H2 on top of graphene and graphite with up to five {Ny) 
vacancies within a total of 120 possible sites of our sim- 
ulation cell. The results, shown in Fig. 3, show a linear 
increase of the superfluid fraction with for both he- 
lium and hydrogen. The increase of superfluidity with 
is significantly larger for ^He; for instance, in the = 5 
case Ps/ p amounts to 0.14 and 0.047 for "^He and H2, re- 
spectively. The effect of the substrate on the superfluid 
fraction is smaller in helium than in hydrogen but in both 
cases Ps/ p is always larger in graphite than in graphene. 
The difference is very small in ^He but significant in H2. 
A possible explanation of the increase of superfluidity 



in graphite relies on the slightly reduced motion in the 
transverse direction due to stronger attraction with re- 
spect to graphene. This increase in confinement makes 
more effective the motion of the center of mass of the 
system in the x — y plane of the simulation box, where 
large particle-permutational rings can be created, and 
consequently Ps/ p increases ([5]). 

We have verified that both in hydrogen and helium the 
periodic spatial order characteristic of the solid phase is 
conserved up to 5 vacancies in 120 possible sites. This 
is shown in Fig. 4, where characteristic snapshots of 
DMC simulations are shown for the perfect and defec- 
tive = 5 film phases (''He). These snapshots corre- 
spond to the probability density that in DMC is repre- 
sented by a collection of walkers, each one of 'iN coor- 
dinates, that evolves in imaginary time according to the 
Schrodinger equation. As one can see, even in the case 
of the perfect crystal there is a finite probability of vis- 
iting the intersite space, which can be interpreted as the 
finite exchange probability leading to non-zero superfluid 
signal (a quantitative estimation of probability exchange 
would require of specific methods beyond the scope of 
the present work^i). When vacancies are present in the 
system the paths connecting different sites become more 
populated, i.e., the superfiuidity increases, but solid or- 
der remains preserved. It is worth noticing that in the 
snapshot of the layer with vacancies one can not allo- 
cate the vacancies implying that they have become in- 
distinguishable. Same snapshots for H2 show that in the 
perfect crystal the intersite occupation is zero (zero su- 
perfluidity) and that the vacancies are identified more 
easily than in helium. 

Summarizing, using the DMC method we have studied 
the supersolidity of the first layer of ''He and H2 adsorbed 
on graphene and graphite in the limit of zero tempera- 
ture. The *He -\/3 x -\/3 commensurate phase shows a 
small but finite superfiuid signal (0.67 %) whereas the 
H2 one does not within our resolution limit (1 x 10~^). 
When vacancies are present in the system, the superfluid 
fraction increases with the concentration of defects; this 
effect is larger in helium where we have obtained values 
as large as ps/ p = 14:%. As the presence of point defects 
in quantum layers is plausible due to imperfections in 
the adsorbent surfaces, further experiments on quasi-2D 
systems can lead to the emergence of new and tunable 
supersolid scenarios i^i^ 
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